D=[-59 -28 0 31 61 92 122 153 184 214 245 275];
t=[9 10.5 12 13.5 15];

for i=1:length(D)
    for j=1:length(t)
        alpha_s(i,j)=alp_s(D(i),t(j));%太阳高度角
        gama_s(i,j)=ga_s(D(i),t(j));%太阳方位角
    end
end